Topological Hamiltonian as an Exact Tool for Topological Invariants 
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We propose the concept "topological Hamiltonian" for various topological insulators in interacting systems. 
Although a bad tool for calculating energy spectra, "topological Hamiltonian" is an exact tool in producing 
topological invariants. We also show that, unlike the topological Hamiltonian, the natural "effective Hamilto- 
nian" for quasiparticles can produce wrong topological invariants, despite its ability to give us excellent energy 
spectra. We present a physical explanation for this counterintuitive phenomenon from the surface state picture, 
more precisely, we show that the topological Hamiltonian determines the existence or absence of gapless sur- 
face modes. This confirms the topological Hamiltonian instead of the mean field effective Hamiltonian as an 
exact tool for calculating topological invariants in interacting systems. This result has significant implications 
for analytical and numerical calculations of topological invariants in interacting systems, in particular, it shows 
that the conventional band structure calculations of topological invariant can fail. 

PACS numbers: 73.43.-f,71.70.Ej,75.70.Tj 
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Introduction. Topological insulators are new states of mat- 
ter whose bulk is gapped just like ordinary insulators, but 
the surface is robustly gaplessHJ-Q]. In the noninteracting 
limit, the gapless surface states are protected by bulk topo- 
logical invariants B5l 4l4ll which are defined in terms of the 
noninteracting Bloch band states. Such topological invari- 
ants are inapplicable to systems with electron-electron in- 
teraction, though many numerical calculations resort to var- 
ious versions of "effective non-interacting Hamiltonian". It 
is worth emphasizing that there is no general justification for 
these extensions of noninteracting topological invariants to in- 
teracting cases, though they may have good chance to pro- 
duce correct result if the interaction is weak. More recently, 
topological insulators with st rong electron-electron interac- 
tion attracts much interest 115143 ill , it is therefore highly de- 
sirable to develop topological invariants for interacting insu- 
lators. Topological field theory description lll ill was devel- 
oped for such interacting topological insulators, however, a 
major question remained to be answered was to obtain ex- 
plicit formulas for the quantized topological responses. It 
was proposed in Ref. B32ll that a Z2 topological invariants 
can be defined in terms of Green's function for interact- 
ing topological insulators in three spatial dimensions (3D). 
There are numerous recent works making wide use of Green's 
function B33l438ll . More recently, exact topological invariants 
in terms of zero frequency Green's function were systemat- 
ically proposedl 391 14011 for various topological insulators in- 
cluding quantum Hall insulators and time reversal invariant 
topological insulators in 2D and 3D, which found interesting 
applications in various contextl 36l 41-451. 



The applicability of the topological invariants expressed 
in terms of zero frequency Green's function is beyond the 
"renormalized bands" (or quasiparticle) picture. Generally 
the quasiparticle picture is not applicable to the calculation 
of topological invariants, because the concept of "quasiparti- 
cle" itself can break down in a strongly correlated system. In 
the present paper, let us focus on systems in which the quasi- 



particle picture works as a good approximation, then one may 
naively expect that an effective Hamiltonian describing quasi- 
particles will produce the correct topological invariants. A 
qusiparticle with an approximate energy u> and momentum k 
has a self energy £(«, k), therefore, if we consider a "flat- 
band" insulator with all filled bands having the same renor- 
malized energy (or quasiparticle spectrum) «o, then we may 
guess that S(a>o,k) will be important for obtaining topologi- 
cal invariants. However, this expectation is incorrect, instead, 
only X(0, k) appears in the exact topological invariants 14011 . It 
is our purpose of the present paper to understand this counter- 
intuitive result. 

The central observation of the present paper is that the zero 
frequency Green's function tells us about the existence or ab- 
sence of the gapless surface states, while the more conven- 
tional "effective Hamiltonian", though an accurate tool for 
bulk energy spectra, fails to do this. We will give an explicit 
example for this surprising phenomenon in the main text. This 
physical observation from surface states confirms the zero fre- 
quency Green's function approach as a crucial tool in topolog- 
ical invariants calculation for interacting systems, which will 
be of great use in searching for new materials for topolog- 
ical insulators in systems with strong electron-electron inter- 
action. We coined the concept "topological Hamiltonian" (see 
the main text) for inverse Green's function at zero frequency , 
from which various topological invariants can be readily cal- 
culated. 

The two main results of this paper can be summarized as 
follows. First, topological Hamiltonian give us bad bulk en- 
ergy spectra (thus it is far from being a good effective Hamil- 
tonian for the bulk), but it gives us exact topological invari- 
ants; second, and more notably, topological Hamiltonian de- 
termines the gapless surface states. 

To conclude this introduction, we emphasize that all sys- 
tems studied in the present paper are gapped in the bulk, 
though the Green's function method is also fruitful in gap- 
less systems[46, 47], in which momentum space defects can 
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be topologically protected. 

Topological Hamiltonian and Effective Hamiltonian. Let 
us fix the notation first. The many body Hamiltonian is 
given as H — Ho + H\, where Ho = 2k<V 1 o(k)ck = 
Tik,afi c l a [ho(^)]a/3Cki3 is the free quadratic Hamiltonian(a,jS 
run through all degree of freedom other than momenta, and 
the chemical potential has been absorbed into the definition 
of Ho) , and Hi term describes the electron-electron interac- 
tion. Various (Matsubara, Retarded, etc) Green's functions 
are defined in the frequency-momentum space in the usual 
way Hit]. For instance, the free fermion system has Matsub- 
ara Green's function G{icD, k) = l/[ia> — ho(k)]. For inter- 
acting systems, Green's functions are still determined by the 
many body Hamiltonian, though their explicit forms can be 
very complicated. 

Recently, in Ref. ll39l l40ll various exact topological invari- 
ants expressed in terms of zero frequency Green's function 
were systematically defined and derived. The idea is outlined 
as follows. The zero frequency Green's function G(iu> = 0, k) 
is a Hermitian matrix, therefore, in the eigenvalue equation 



G _l (0, k)|a(0, k)) = fi a (0, k)|a(0, fc)> 



(I) 



all /u a (0, k) are real numbers. An eigenvector |ff(0, k)} with 
yu„(0, k) > is called an "R-zero", and the space spanned by 
R-zeros are called "R-space". We can define Berry gauge field 
and Berry curvature in the R-space, and then define topologi- 
cal invariants just in the same way as the noninteracting cases. 
For example, for two-dimensional topological insulators with 
broken time reversal symmetry (namely the quantum Hall in- 
sulator), we have the "generalized Chern number"[40] 



Ci 



2nj 



d 2 kT, 



(2) 



where Tij - diJ[j - djMi, and Mi - -i ^(kald^lkar). Here 
|ka) runs through an orthonormal basis of the R-space. This 
simple generalized Chern number was shown j40ll to be as 
powerful as the more complicated topological invariants de- 
fined using frequency integral||32L |4q, |49|] . Systematic gener- 
alizations of this zero frequency Green's function formula to 
topological insulators and superconductors in other symmetry 
classes can be found in Ref. lHOi 14211 . and certain special case 
was studied in Ref.[50[ 51 ] from different motivation. 

Since G l (Q, k) plays an crucial role in our approach, let us 
define the "topological Hamiltonian" as 



ft topo (k) = -G- 1 (0,k) 



which can also be written as 



Wk)^o(k) + S(0,k) 



(3) 



(4) 



following the Schwinger-Dyson equation. With this defini- 
tion, the generalized Chern number given in Eq.® is equiva- 
lent to the Chern number of an noninteracting system with free 
Hamiltonian ho(k) = /z topo (k). Other topological invariants 



such as Z2 invariants can also be calculated from the topolog- 
ical Hamiltonian in similar way|40]. An naive understanding 
of /itopo is to think it as an effective Hamiltonian for quasipar- 
ticles, namely that h topo is a formula for "renormalized energy 
band". As we will see, this understanding is incorrect. In fact, 
^topo is far away from being an effective Hamiltonian, which 
is expected to properly describe the quasiparticle. One of the 
failures of h topo to be an effective Hamiltonian is that its eigen- 
values can be far away from the true energy spectra, as we will 
show below. The powerful aspect of topological Hamiltonian 
approach is its ability to produce topological invariants rather 
than energy spectra. 

What is the natural effective Hamiltonian producing ac- 
curate energy spectra? Let us start from the self-consistent 
equations ll52ll53ll for quasiparticle spectra, which can be writ- 
ten in a simplified notation as 



[h (k) + S(w a ,k)]|K ff (k)> = w (k)|ii (k)> 
which can be more compactly written as 
G~\cj a ,k)\u a (k)) = 



(5) 



(6) 



Eq.© is a self-consistent equation with the self energy 2 de- 
pending on the spectra w„. It is worth to emphasize that 
it is 2(Wq., k) instead of E(0, k) that appears in Eq.©. In 
fact, this simple observation was the initial motivation of the 
present paper. In a general interacting system, E(0, k) can be 
far away from li(o) a , k), therefore the topological Hamiltonian 
approach using X(0, k) is a bad tool for energy spectra, though 
it is an powerful tool for topological invariants. A natural 
question is whether we can define physically meaningful topo- 
logical invariants in terms of S(w ff ,k) or G _1 (w ff ,k). As we 
will see in the following, this effective Hamiltonian approach 
turns out to be less fruitful than the topological Hamiltonian 
approach. 

Let us proceed with this 2,(cj a , k) approach to topological 
invariants. It is natural to define a frequency dependent "ef- 
fective Hamiltonian" 



/i eff («,k) = /jo(k) + X(«, k) 



(7) 



Furthermore, it seems also natural to define Berry connection 
in terms of \u a (k)) instead of |a(0, k)) given in Eq.([T]l, since 
|w a (k)) more faithfully represents the quasiparticle picture in 
the interacting system. Now there is a difficulty in doing so, 
namely that k) is generally not a Hermitian matrix, but 
we will show that even if this difficulty can be circumvented, 
the topological invariant defined in this approach is still in- 
correct. Let us write £(o>, k) = 2,\(u>, k) + H^(a),1t), with 
both Ei and £2 Hermitian. The assumption that quasiparti- 
cle picture is a good approximation for our system is equiv- 
alent to saying roughly that £2 can be ignored compared to 
Si, since a non-vanishing S2 gives an imaginary part to the 
quasiparticles and implies finite lifetimes. Good quasiparticle 
approximation implies that the the inverse lifetime should be 
significantly smaller than the characteristic energy scale of the 
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quasiparticle. Therefore, let us replace 2,(a>, k) by ~Li(a>, k) in 
Eq.©, namely 



[h (k) + £i(w a ,k)]K(k)> = w„(k)|M„(k)> 



(8) 



Here we encounter another difficulty, namely that \u a {k)) are 
generally not orthogonal. This is a relatively minor prob- 
lem since we can obtain a orthonormal basis by the Gram- 
Schmidt process. For 2D quantum Hall systems, let us define 
the "quasiparticle Chern number" as 



2nJ 



d 2 kT„ 



(9) 



where Ttj = d^j - djM,, and ^ = -i2 a (S Q ;(k)|fl^|M Q .(k)), 
in which |5 tt (k)) is an orthonormal basis in the space spanned 
by those |w tt (k)) with negative co a (k). For the simple model 
to be studied in the next section, \u a (k)) themselves will be 
sufficient. 

In the remaining parts of this paper, we will show by ex- 
plicit example that Eq.® differs from the correct topological 
invariants given by Eq.®, therefore it must be incorrect. We 
then explain the reason why Eq.© instead of Eq.© is the 
correct topological invariants, from the surface state picture. 

Effective Hamiltonian can produce wrong topological in- 
variant. In this section we will give an explicit example in 
which the effective Hamiltonian approach to topological in- 
variants is incorrect, for instance, Eq.© can fail for quantum 
Hall system even if the quasiparticle picture is good. Let us 
begin with a two-dimensional (2D) model of quantum anoma- 
lous Hall effect ll54l 15511 with electron-electron interaction. Let 
us assume the free Hamiltonian part Ho describes a two band 
model[55], with 

ho(k) = sin k x cr x + sin k y <T y + [mo + (2 - cos k x - cos & v )]oty.O) 

Instead of giving explicit form for electron-electron interac- 
tion terms H\, in the following we will choose an ansatz for 
self energy l,{a>, k), which is sufficient for our purpose. 

Let us make a simple ansatz that k) ss £ z (w, k)cr, near 
k = 0, in which case the main idea can be most readily ap- 
preciated. Near k = we have h e ff (at, k) m k x cr x + kyCr y + 
[mo + k = Q)]cr z . Additionally, we assume that the inter- 
action preserves particle-hole symmetry of the free fermion 
Hamiltonian Ho, implying that 0) = Y, z (-u), 0). From 
Eq. ©(taking k = 0) we can obtain the self-consistent equa- 
tion 



ffief = mo + S,(m eff , k = 0) 



(11) 



for the two eigenvalues +m e g of /i e ff(w, k = 0). This effective 
mass will be crucial in calculating C\. On the other hand, the 
topological mass is defined as 



the topological Hamiltonian /i topo is given as 

Ci = f efkn ■ (d k n X d k A) (13) 
4n J 

where the integrating range is the Brillouin zone, n = 
(n x , n y , n z ) is extracted from h lopo = cr ■ n, and n = n/|n|. As a 
concrete example, let us pick mo = -0.2, £,(0, k = 0) = 0.1, 
then m lopo = —0.1, from which it follows that C\ = -1. If 
£,(0.1, k = 0) = 0.3 (as an example), then m s s = 0.1 is 
self-consistently obtained from Eq.dTTI). It follows that C\ 
determined from Eq.® is C\ = 0. We knew that Eq.(f2| 
produces the correct topological invariants C;|40], therefore, 
Eq.© fails in this case because 



Ci *C X 



(14) 



"Jtopo = m Q + 1,(0, k = 0) 



(12) 



Topological Hamiltonian determines topological surface 
states. In the previous sections we have shown that the mean 
field calculation with the natural effective Hamiltonian /i e ff 
produces wrong topological invariants, on the other hand, the 
"topological Hamiltonian" h topo (k) = -G~'(0, k) produces 
the correct topological invariants. This may seem mysterious 
since h topo is a bad tool for quasiparticle spectra. Why should 
we expect /i topo instead of heg to produce the correct topolog- 
ical invariants? The answer, as we will show, is the fact that 
the gapless surface states are determined by h topo . The key 
idea is that the zero energy surface states, if exist, will "feel" 
the self energy l,(ai = 0) because they have zero energy. In 
sharp contrast to /i topo , h e ff is not directly related to the exis- 
tence or absence of surface modes, despite of the fact that they 
they tell us the energy spectra of the bulk . As a well known 
fact, the existence of topologically protected surface states is 
a holographical manifestation of the bulk topology, therefore, 
the connection between /z topo and surface modes firmly estab- 
lishes the h topo approach to topological invariants. 

Now let us flesh out the above idea in the simple two- 
band model studied in the previous section. To this end, let 
us first assume that there are zero energy surface states, and 
then see what it implies. Consider a 2D system with coor- 
dinate (jc,y)[see Fig.©]. What we shall show is that at the 
interface between two bulk states with different C\ (calcu- 
lated from topological Hamiltonian), there is gapless surface 
state. Let us suppose that the many-body Hamiltonian H is 
a function of x, namely that certain parameters (e.g. inter- 
action strength) of H vary as functions of x, then the topo- 
logical Hamiltonian h topo is also a function of x. Assuming 
the form of /i top o given in the previous section, we now have 
h topo (x) a —icr x d x +cr y k y +m topo (x)cr z since /i topo is now a func- 
tion of x, and the translational symmetry along x-direction is 
broken. Because it will not affect our conclusion, we take 
l^topoMI << 1 to simplify our analysis. Let us suppose 

witopoW M (x» 0); m topo (x) -> -M (x « 0) (15) 



Neark = 0, we have h topo (k) m k x a ' x +k y o- y +mt apo o~ z . Follow- 
ing Eq.©, the correct topological invariants calculated form 



and w topo smoothly interpolate between this two limits near 
jc = 0[seeFig.CQ>)]. 
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Zero Mode 




FIG. 1: The gapless mode on the edge, (a) The physical setting and 
the zero mode, (b) m topB as a function of x. 



Due to the translational symmetry in y direction, k y is a 
good quantum number, thus we can study the k y = branch 
first. As we have emphasized, the topological Hamiltonian 
/i topo (k) = -G-\0,k) = ho(k) + S(0,k) is a bad effective 
Hamiltonian for the bulk because the true energy spectra are 
nonzero in the bulk; on the other hand, /i topo is the tailor-made 
effective Hamiltonian for the zero mode living at the inter- 
face. Because the zero mode has exact zero energy, the self 
energy is given exactly by E(w = 0). In other word, locally the 
zero mode is described by an "effective Hamiltonian" /i t0 po(*)- 
With this key observation, it is now most natural to study the 
zero mode solution as the eigenvalue problem of /i t0 po- Ex- 
plicitly, we have 



h 



topo 



-icr x d x + 0- z m topo O) 



(16) 



where we have set k y = as mentioned before. We have 
omitted terms in higher order of k since we are concerned 
with low energy modes. If M > [M is defined in Eq.(fT3bl. 
the solution of zero mode from Eq.dToTi is readily obtained as 



\i]f{xj) 



zero mode 



dx'm topo (x')] (17) 



The assumption of existence of this zero mode requires that 
(x) < when x — > -co and m tnri Jx) > when x — > oo 



Hopo 



(otherwise the wavefunction diverges)[56, 57], which is satis- 
fied by this solution. If M < 0, then there is also a zero mode 
solution with (l,-i) T replacing the (l,z) T factor in Eq.([T7li. 
Note that the existence of zero mode is fully determined by 
the asymptotic behaviors of h topo (x) as \x\ » 0, therefore it 
is determined by the bulk topology. From the above explicit 
calculation at the interface we have seen that when m topo (x) 
changes sign at the interface near x — 0, there are exact zero 
modes. The above calculation looks in parallel with the non- 
interacting cases in which we start form /iq instead of /i topo , 



however, the crucial difference is that in interacting cases we 
have a frequency-dependent self energy, and "£(u> = 0) appears 
in a self-consistent manner. 

From the bulk calculation given in the previous section, we 
have Ci(M > 0) = and Ci(M < 0) = -1 (remember that 
\M\ « 1), thus the topological invariants at the two sides of 
interface are different. To summarize the above calculations, 
we have shown in an explicit example that when the topolog- 
ical invariants calculated from h iopo are different at the two 
sides of an interface, there are gapless surface states. In ex- 
periment, one side of the interface is often the vacuum, which 
can be regarded as a topologically trivial insulator, therefore, 
"different topological invariants" at the two sides of interface 
means that the insulator has nontrivial topological invariant 
(Ci + in this quantum Hall example). The fact that topolog- 
ical Hamiltonian /i topo (k) determines the existence or absence 
of gapless surface mode is among the central result of this 
paper. It is equally notable that /i e fF plays no role in the solu- 
tion of zero energy surface mode. Shortly speaking, we have 
"hes for the bulk states and /i t0 po for the surface states", though 
both /left and /i e ft are defined in the bulk. This is the deep rea- 
son why /i topo instead of h e s determines the bulk topological 
invariants. 

Conclusions. In conclusion, we have shown that the "topo- 
logical Hamiltonian" /i topo , defined in terms of zero frequency 
Green's function, faithfully produces the topological invari- 
ants for interacting topological insulators, while the natural 
qausiparticle effective Hamiltonian fails to do this, though 
it gives much better energy spectra. The deep reason, as we 
have shown, is that /z topo directly tells us the information of 
gapless surface states, while fails in this respect. Although 
we only write down explicit formulas for 2D quantum Hall in- 
sulators, the same analysis can also be applied to topological 
insulators in other symmetry classes. It is also worth noting 
that the purpose of this paper is to propose and understand 
the topological Hamiltonian approach as an exact tool, while 
more realistic models in which this approach will find wide 
applications remain to be studied. 

Let us conclude with the remark that in band structure cal- 
culations we are generally more concerned about quasiparticle 
spectra, namely that various versions of the effective Hamil- 
tonian approach are preferred, but as we have shown in this 
paper, such calculation can fail in topological invariant calcu- 
lations. The topological Hamiltonian approach, on the other 
hand, faithfully produces correct topological invariants, thus 
it is an crucial tool in searching for topological insulators with 
strong electron-electron interaction. 
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